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Abstract 

We study analytically the effect of temperature cyclings in mean-field spin- 
glasses. In accordance with real experiments, we obtain a strong reinitialization 
of the dynamics on decreasing the temperature combined with memory effects when 
the original high temperature is restored. The same calculation applied to mean- 
field models of structural glasses shows no such reinitialization, again in accordance 
with experiments. In this context, we derive some relations between experimentally 
accessible quantities and propose new experimental protocols. Finally, we briefly 
discuss the effect of field cyclings during isothermal aging. 



LPT-ENS/9842; LPTHE/9853. 



Glasses are characterized by having extremely slow relaxations and by the strong 
dependence of their behaviour upon the ( "waiting" ) time elapsed since their prepa- 
ration. The latter property is usually called physical aging. 

A means to study the dynamics in the glassy phase in more detail consists in 
following the evolution of the sample under a complicated temperature history. The 
protocols that have been more commonly used include temperature cyclings within 
the low temperature phase. 

The results for different types of glasses are quite different .0~i Spin-glasses show 
the puzzling phenomenon of reinitialization of aging following a decrease in temper- 
ature, combined with the recall of the situation attained before the negative jump 
when the original high temperature is restored.!!! Remarkably, when similar proto- 
cols were applied to structural glasses, e.g. in dielectric constant measurements of 
glycerol by Leheny and Nagel,! no substantial reinitialization was observed.0 

This difference in the effect of temperature changes on spin and structural glasses 
is a fact that any generic theory of glasses is expected to explain. 

Different groups interpreted the behaviour of spin-glasses under temperature cy- 
clings during aging as evidence for both the droplet&i and the hierarchical! pictures 
of the dynamics (the former with some extra refinement! with respect to the original 
versions of the eighties^). 

The hierarchical dynamic picturelll'! is a heuristic way to think about the results 
from positive and negative cyclings inspired in the organization of equilibrium states 
in the Parisi solution of mean-field spin-glasses, such as the Sherrington-Kirkpatrick 
model. It is assumed that spin-glasses have a large number of metastable states 
that are organized in a hierarchical fashion just like the equilibrium states. It is then 
proposed that the system is composed of (independent) subsystems whose dynamics 
is given by the wandering in such a landscape. An average over subsystems has to 
be invoked in order to obtain smooth results as observed in experiments. 

A concrete realisation of a hierarchical dynamic system can be made with the trap 
models.EUll These models have been solved analytically in isothermal conditions.@ 
Though a full analytic description of their dynamics in the presence of tempera- 
ture cyclings is not available yet, a careful discussion of their effects yields very 
encouraging results.El 

Surprisingly enough, the main features of the cycling experiments have never 
been derived analytically from microscopic models, while the numerical evidence&izl 
is inconclusive. In this paper we shall show analytically how these effects arise in 
mean-field models of spin glasses, and why they are absent in mean field models 
of structural glasses. One of the questions that will receive a clear answer is why 
the effects should be hardly observable at very short times, such as are inevitably 
involved in simulations. 

We shall consider in detail the particular class of temperature cycling experi- 
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FIG. 1. The TRM curves after a long negative cycle sketched in the inset. The tem- 
perature at which the relaxation is measured is T = 12K = 0.72 T g . The bold lines 
are reference curves corresponding to waiting at a constant temperature of 12K during 
t w = 10, 30, 100, 1000 min. The thin lines are obtained after a temperature cycling of 
1000 min and temperature jumps AT = 0.1, 0.3, 0.5, IK. The data is courtesy of E. 
Vincent, see Ref. [2]. 

ments in which the thermoremanent magnetization (TRM) is measured. Similar 
conclusions havebeen extracted from the out of phase susceptibility (%") data at 
fixed frequency.@c3 

There is however a slight difference between TRM and ac x" measurements. 
In the former the TRM after the cycling is recorded and, since measurements are 
directly performed in the time domain, one has access to very large time scales after 
t w . In the latter case the x" is measured during and after the cycling This allows 
us, for example, to clearly see the large reinitialization of the dynamics provoked by 
the negative jump.ll The price to pay is that in ac measurements the frequencies are 
necessarily small compared to to the inverse of the measuring time. One then has 
access to relatively small time-differences. 

In the TRM experiments of Refregier et al.i the system is quenched to a subcrit- 
ical temperature under a small field that is used as a probe, for which the linearity 
in the response is checked within the same experiment. 

In the negative temperature cycling experiment (see the inset of Fig. [I]), the 
system is quenched at t = to a temperature T + . At a time t\ the temperature is 
dropped to T~ and at a later time t<i the temperature T + is restored. The fraction 
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FIG. 2. The TRM curves after a short positive cycle sketched in the inset. The 
relaxation is obtained at T = 12K = 0.72T g . The bold lines are reference curves obtained 
after isothermal waiting at 12K for t w = 30, 1000 min. The thin lines are the decay of 
the TRM after a cycle of 5 min and AT = 0.25, 0.5, 0.75, 1, 1.5, 2.5 K. Data is courtesy 
of E. Vincent, see Ref. [2]. 



of time spent at T~ is rather large. 

The resulting state of the sample after the temperature cycling is investigated by 
cutting off the field at a time t w soon after t 2 and by recording the subsequent decay 
of the TRM. The results for different AT are shown in Fig. |l|: The TRM curves for 
each AT can be superposed to the TRM curves obtained at constant T + but with 
an effective waiting time t\ + t w — t 2 < £ FFF (AT, t\, t%, t w ) < t w . Note that the first 
inequality implies that the system remembers the evolution performed at the higher 
temperature T + while the time spent at the lower temperature T~ is partially (or 
even totally) erased. 

Once the negative cycle is over, its main effect is to slow down the aging process. 
This result is very intuitive and will hold for almost any system with slow dynamics 
activated by thermal noise. 

The real surprise appears when a cycle of positive temperature (inset of Fig. ||]) 
is applied to spin glasses. Here the procedure is the inverse: the system is quenched 
to a temperature T~ up to a long time t\. This is followed by a short period at 
temperature T + from t\ to t 2 , at which T~ is restored. As before, the field is cut 
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off at a later time t w and the subsequent decay of the TRM is recorded. 

The result is shown in Fig. |2|. The higher the upward pulse in temperature, the 
generally younger the system seems but, unlike in the negative cycling case, the 
effect cannot be described with an effective waiting time. 

Note that in the simplest cases of aging through activated processes (as for 
example the coarsening of the random field Ising model) the effect would be the 
reverse: the temperature pulse would quicken the activated processes and help aging. 
Refreshment would arise only if the temperature pulse is high enough to take the 
system above the transition. 

Indeed, the solution for mean-field spin-glasses discussed below involves an addi- 
tive separation of the TRM curve M TRM = M REP + M F . In the positive temperature 
cycling the "slow" component M REP is not affected by the pulse while the "fast" 
component M F is taken to its effective critical temperature and is thus completely 
reinitialised. The amplitude of M F (M REP ) is larger (smaller) for higher (lower) tem- 
perature pulses. In the negative temperature cycling one is essentially probing only 
the slow component M REP that is slowed down by the effect of the lower temperature 
excursion. 

We shall show that the difference in the effects of positive and negative tem- 
perature changes is present in the mean-field version of spin-glass models while it is 
absent in the models thought to be a mean-field caricature of structural glasses. The 
relevant difference between these models resides in their dynamic behaviour below 
the transition. The former decay in infinitely many time-scales while the latter do 
in only two. (For an unambiguous definition of time-scales in aging problems see 
Ref. fl6|1 .) The explanation we give here is based on this difference. We choose to 



discuss in detail the magnetization behaviour although the susceptibility data, in 
particular the reinitialization after a negative jump, can also be understood within 
this framework. 

Another common way to test the dynamics in the glassy phase is to apply mag- 
netic field steps during aging at constant temperature. In these experiments no 
important difference is obtained between switching on and switching off the dc field. 
We shall briefly discuss this result within the same analytic framework. 

The organisation of the paper is the following. In Section [I] the main features 
of the analytic approach are described. In Section U the temperature cycling ex- 
periments are explained within this analytic approach. The results of field cycling 
experiments are confronted to this approach in Section |TJ. Section [TV] is devoted to 
some new experimental proposals and Section [V] to the conclusions. 
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I. THEORETICAL APPROACH 



In order to understand the effect of temperature cyclings on the relaxation of 
the TRM and out of phase susceptibility of mean-field spin-glass models we have to 
understand the time- dependence of the auto-correlation and response functions: 



1 N I N 



5( Si (t)) 



(1.1) 

h=0 



For mean-field models, one obtains a set of coupled equations that entirely determine 
the dynamics of these two-point functions. The thermoremanent magnetization 
M TRM (i, t w ) at time t, after cutting off a small field at time t w , is then expressed in 
terms of the integrated susceptibility x(t,t w ): 

X (t,t w ) = f drR(t,r) , M TRM (t,t w ) = h[ X (t, 0) - X (t, *«,)] (1.2) 
with h the strength of the small dc field applied. 



A. The model 



For definiteness let us consider the toy model consisting of N continuous spins 
with a spherical constraint Yli s t{t) = N and a random energy V(s) correlated 
as V(s)V(s') = Nu(s ■ s' /N). The statics^ as well as the constant temperature 
dynamics of this model have been solved in all detail (see Refs. [19|~21[|). Two types 



of models with potential correlations such that \/Jv"{C) is, for all < C < 1, 
concave (convex) yield very different statics@ (one step replica symmetry breaking 
versus full replica symmetry breaking) and dynamicstj (first order versus second 
order dynamical transitions) and correspond to mean-field versions of structural 
and spin glasses respectively. Examples that have been extensively studied in the 
literature are the u p = 3 spherical model" for glassesjl! with u(C) = 1/2C 3 (and 
(v"{C))~ l l 2 concave) and the u p = 2 plus p = 4 model" for spin-glasses,il with 
u(C) = l/2(c 1 C 2 + c 2 C 4 ) and c\ and c 2 two constants such that (v"(C))~ 1/2 be 
convex. 

The exact equations of motion for C and R at times t > t' are!l] 
dC(t, t') 



dt 
dR(t, t') 



-z(t)C(t,t')+ f drD[C{t,T)]R{t',r)+ f dTE[C{t,r)}C(r,t') , (1.3) 
Jo Jo 

-z{t) R(t, f) + f dr E[C(t, r)] R(t, f) , (1.4) 
Jt 1 

z{t) = T(t) + /* dr (D[C(t, r)]R(t, r) + E[C(t, r)]C(t, r)) , (1.5) 
Jo 
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with £(t,r) = D'[C{t,r)}R{t,T), D'[C] = d c D[C] and D[C] = v'{C). The function 
z(t) is a Lagrange multiplier that enforces the spherical constraint. We shall con- 
centrate on the spin glass-like case that corresponds to a concave random energy 
correlation. We shall briefly discuss the structural glass-like case at the end of this 
Section and explain why it does not show large asymmetry effects. Let us recall the 
constant temperature solution for the mean-field spin-glass models. 

One of the salient features of the relaxation of mean-field spin-glasses below the 
dynamic critical temperature Td is the presence of infinitely many (two-time) scales 
organised in a hierarchical way. In the low temperature phase the two-point corre- 
lation and response depend on the two times involved. The form of the relaxation 
is usually analyzed by looking at these functions at fixed (but large) waiting time in 
terms of the time difference r = t — t w . The correlation and response functions have 
a first fast stationary relaxation; for instance, the correlation rapidly decays from 
1 to qEA, the Edwards-Anderson parameter. This time-scale is usually called the 
FDT regime, for reasons that will become clear below. For longer time-differences 
the relaxation continues at a waiting-time dependent speed. Furthermore, if one 
imagines the subsequent decay of the correlation as taking place in infinitesimal 
steps, each step takes much longer than the previous one. Indeed, in the limit of 
large waiting time each infinitesimal step implies a period of time that is infinitely 
longer than the previous one and these time-scales get completely separated. The 
same separation of time scales characterizes the decay of the response. 

The sharp separation of time-scales allows us to split the correlations in a fast 
part C F (t,t'), going from one at equal times to go a t very distant times, and a slow 
part C qo (t,t'), going from C qo (t,t) = g and tending to zero at even more distant 
times. The point at which we split the correlation is chosen arbitrarily provided it 
satisfies go < Qea- Correspondingly we separate the response in a fast and a slow 
part. We then have 



Since C qo (R qo ) is infinitely slower than C F (R F ) their time evolution can be charac- 
terized as follows: for all times such that C F changes C qo is just constant and equal 
to go- Instead, in the time-regime in which C qo varies C F has achieved its asymptotic 
value g and does not further evolve. 

Note that the time-scale separation is achieved only in the large waiting-time 
limit. This is the crucial ingredient for the argument we shall use to show that 
Eqs. ( |1 . 3| )- (p~T5|) capture the phenomenology of temperature cyclings. We believe that 
the effects disussed in the Introduction have not been clearly observed numerically 
because the times explored were inevitably very short and the time-scales could not 



C(t,t')=C qo (t,t')+C F (t,t')-q , 
R{t,t') = R qo {t,t') + R F (t,t') . 



(1.6) 
(1.7) 



be sufficiently separated.0~§! 
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FIG. 3. Sketch of the % vs C plot. The two straight lines have slopes — 1/T + and 
— 1/T~ and correspond to the FDT result. The curve is temperature independent and it 
is given by (v" {C))~ l l 2 . The bold segment is the part erased (created) when decreasing 
(increasing) the temperature. See text for more details. 

As we shall show, the change in time-dependence of the fast and slow parts of 
C and R under a temperature jump are very different. The slow parts C qo and 
R qo are modified through a smooth time-reparametrization. The fast parts C F and 
R F behave as the correlation and response of an effective model with a critical 
temperature precisely equal to T + . The temperature jumps have a strong effect on 
the fast parts and this effect is very different depending on the sign of the change. 

The other salient feature of the low-temperature dynamics of mean-field spin- 
glasses is the violation of the fluctuation-dissipation theorem (FDT) and, most im- 
portantly, the generalized form that the relation between correlation and response 
takes. 

A very useful way to quantify the fluctuation-dissipation relation that holds out 
of equilibrium is given by the relation between the integrated response x(Mu>) and 
the correlation C(t,t w )M At fixed and large t w , one constructs a plot of vs 
C(t,t w ) using t as a parameter. This is shown in Fig. |3| for two temperatures T + 
and T~ . For each temperature the x vs - C curve consists of two parts: 

• A straight line of slope minus the inverse temperature. This corresponds to 
the fast time regime where the FDT holds and C decays from 1 to qEA- We 
call this time-regime the FDT regime. 



• A curve given by 



f(C) 



1 



(1.8) 



v"{C) 



This corresponds to the slower time-regimes where the FDT is violated and C 
decays from q^A to in a waiting-time dependent manner. 

For temperature T + the x vs C pl°t follows the straight line of gradient — 1/T + 
from (1,0) to (qo,Xo) an d then the curve f{C) up to (0, %°°). For temperature 
it is given by the straight line of gradient —1/T~ from (1,0) to (qi,Xi) an d then 
the curve f(C) up to (0, x°°). The Edwards- Anderson parameters are q\ and q , 
respectively. The two extreme cases are T = and T = T c (the critical temperature). 
The former corresponds to a vertical line starting at (1, 0) that matches the curve 
f{C) at C = 1 and then follows it up to (0, x°°). The x vs C plot for the critical 
temperature is just a straight line linking (1,0) and (0,x°°) with a slope — 1/T C . 

The remarkable fact of this family of models is that the curved segment from 
(q , Xo) to (0, x°°) is the same for the temperatures below T + . This is not a general 
property of mean-field spin-glass models. It holds exactly for the family of models 
here considered but only approximately for the Sherrington-Kirpatrick model. In the 
static replica approach, the corresponding property of the function x(q) is called the 
"Parisi- Toulouse approximation" .E3 The numerical evidence seems to show that this 
approximation works very well in finite dimensional spin-glasses,il at least within 
numerically accessible times. 

The analytic argument we shall develop below can be most easily pictured by 
considering Fig. |3|. The range of correlation and response values in which the sys- 
tem has aging dynamics is given by the curved part of the plot. On changing the 
temperature, the straight line corresponding to the fast relaxation moves clockwise 
and anticlockwise like a windshield- wiper, creating and destroying the bold segment 
(<?cb Xo) - ► (<7i>Xi); thus restarting and erasing the aging scales corresponding to this 
interval. 



In the limit of large waiting times we can separate the different time-regimes as 
follows. The equations for the fast parts of the decay are 



B. Analysis 



dC F (t,f) 
di 



-(z F (t) + z(t))C F (t,t')+ f drD[C F (t,r)} R F (t',r) 



Jo 



+ f dr D'[C F (t,r)} R F {t,T)C F (r,t') + z{t) , 



Jo 



(1.9) 
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^Ml = -( z ^t) + z(t))R F (t,t') + J t t i dr D'[C F (t,r)}R F (t,T)R F (r,t'), (1.10) 

z F (t) = T(t) + f dr(D[C F (t,r)] + D'[C F (t,r)]C F (t,r)) R F (t,r) . (1.11) 
Jo 

The equations for the slow parts of the decay are 

small = - (z F {t) + z(t)~ J^drD'[C F (t,r)}R F (t,T^j C qo (t,t') 

+ X F (t')D[C qo (t,t')]+ f drD[C qo (t,T))R qo (t',r) 

Jo 

+ f drD'[C qo (t,T))R go (t,T)C qo (T,t') , (1.12) 

small = - (z F (t) + z(t) - D'[C qo (t,t')]x F (f) ~ f dr D'[C F (t,r)]R F (t,r)j R qo (t,t') 

+ £drD'[C qo (t,r)}R qo (t,r) R qo (r,t') , (1.13) 

z{t) = fdr (D[C qo (t,r)} + D f [C qo (t,r)]C qo (t,r)) R qo (t,r) . (1.14) 



The input of the slow scale on the fast one is given by the one-time quantity z(t). The 
fast scale influences the slow one through three factors: z F (t) defined in Eq. ( |1 . 1 1| ) , 
X F {t) an d X F (t) defined as 

X F {t)= fdrR F {t,T), 
Jo 

POO 

X F (t) = jf drR F (r,t), (1.15) 

and the last term in the parenthesis of Eqs. ( |1.12| ) and ( |1.13| ). One can eliminate 
the dependence on this last factor by using the relation 

small ~ - (z F (t) + z(t) - J* dr D'[C F (t, r)]R F (t, r)j q + X F (t) D[q ) + z(t) (1.16) 

that follows from Eq. ( |1.12j ) evaluated at 

t'^t C qo (t,t') -> q . (1.17) 

After inserting the relation ( |1.16| ) in Eqs. ( |1.12|) and ( |1.13|) one has the following set 
of slow equations 

sma ll = --(z(t) +x F (t)D[q ])C qo (t,t') + f dr D[C qo (t,r)] R qo (t',r) 
qo Jo 

+ X F (t')D[C qQ (t,t')] + fdr D'[C qo (t,r))R qo (t,r) C qo (r,t') , (1.18) 

Jo 
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small = - L(z(t) + x F (t)D[q ])R qo (t,t') + 5ht* \t')D'[C qo (t,t')]R qo (t,t') 
Qo 

+ f dr D'[C qo (t,r)]R qo (t,r) R qo (r,t') , (1.19) 

together with Eq. ( |1.14| ) for z(t). The effect of the fast scales on the slow one is 
then only given by z F (t), x F if) an d X F (^)- 

At constant temperature the quantities z(t), X F {t) an d X F (t) tend to the following 
limits 

z(t)^z°°, X F (t)^Xo, X F (t)^Xo- (1-20) 

As a consequence of the independence of temperature of x(C) if C < Qea-, see Fig. ||, 
neither z°° nor xo depend on the temperature. 

Consider now the temperature cycling experiment. We consider the effect on the 
slow and fast parts separately and denote the correlation after the cycle 

C CYCLE (t, t') = C« FP (t, t') + C F (t, t') . (1.21) 

The first term is the slow part while the second is the fast part. Similarly, for the 
TRM: 

M CYCLE (t, f) = M™ p (t, t') + M F (t, f) . (1.22) 

Under the assumption that in the presence of the cycling C qo and R qo are still 
slower than C F and R F and are affected at most by time-reparametrisations, an as- 
sumption whose consistency has to be checked, z(t) remains constant and equal to z°° 
throughout the cycle. The function x F (^) stabilises to the temperature-independent 
value Xo soon after each temperature change on a time-scale that is very short com- 
pared to any aging times since it is the relaxation of a one time-quantity.^ Apart 
from a short pulse, proportional to x{t) ~ Xo, the slow equation feels the effect of 
the temperature cycling only through the small matching terms, that we collected on 
the left-hand-side of Eqs. ( |1.12|) and ( |1.13|) under the name of "small". These affect 



only reparametrizations,0 consistently with the initial assumption in the discussion 
Hence, the slow parts of the correlation and response are modified by the cycling 
through a smooth change in the speed of evolution, i.e. a time-reparametrisation. 
This is more precisely stated by studying the change of the slow parts in each time- 
scale. With this purpose, let us denote C^ + (t,t') and R^, + (t,t') the correlation 
and response function in the time-scale labelled a for a system aging at constant 
temperature T + . In the presence of a temperature cycle between ti and t 2 each of 
these is modified according to 

cr CLE (M') = C T a + {h a {t) : h a {t')) , iC CLE (M') = d t rh a (t') R T a + (h a (t),h a (t')) . 
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Though the calculation of h a (t) is beyond our analytical means, we know that it 
has three distinct behaviours depending on the relation between t', t and ti, t 2 . 
Obviously, h(t) — t if t < t%; h a (t) has inflections around t\ and t 2 ; and we expect 
dth a (t) ~ 1 for t ^> t 2 when the svstem forgot the cycle, in accordance with the 
weak long-term memory property.00 Furthermore, we expect that during the period 
t\ < t < t 2 the dynamics slows down, dth a (t) < 1, if we cool down the sample, while 
the dynamics accelerates, dth a (t) > 1, if we heat the sample. In Eq. ( |1.21| ) we have 
generically denoted the reparametrized slow part C™ F . 

The fast equations (|1.9|)-( |i"TT| ) decouple from the slow equation and correspond 
to an effective model in a constant field, since one can check that the (essentially 
constant) term z(t) that appears in the fast equation can be looked upon as the 
effect of an effective magnetic field. In fact, T + is the critical temperature of the 
model in the presence of such a field. 

The temperature of the effective fast model is cycled following any of the two 
experimental protocols between its subcritical temperature T~ and its critical tem- 
perature T + . This is why the effect of the temperature change is stronger on the 
fast part of the correlations and heavily depends on the sign of the change. 

The structural-glass models in which (y"{C))~ l l 2 is concave have a quite simpler 
isothermal dynamic behaviour below its transition. They decay in only two two- 
time scales with a first stationary decay towards q^A and a subsequent slower non- 
stationary decay towards zero. For these models we are not allowed to separate the 
correlation and response in the additive way proposed in Eq. ( |1.7| ) for any q and 
all the analysis presented in this seccion simply does not apply. 



II. COMPARISON WITH THE EXPERIMENTAL RESULTS 

Let us now discuss the experimental TRM curves in the light of these results. 



A. Negative temperature cycling 

The straight line in the \ vs. C curve presented in Fig. ^ moves clockwise and 
then anti-clockwise to its original position. At the end of the cycle the bold part of 
the curve disappeared from the x vs C plot. At t = the fast model is quenched 
down to its critical temperature T + , and rapidly equilibrates. At time t x it is 
further quenched to its subcritical temperature T~ and starts aging with the effective 
waiting time t^ FF = t w —t\. Finally, at time t 2 the fast model is brought again to 
the critical temperature T + , and it rapidly reequilibrates. The fast part has only 
made an excursion into its glassy phase through its quench to T~ but, since it has 
been taken back to its critical temperature T + , it has quickly forgotten it.0 
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One concludes that there is then an effective waiting-time t FFF 

tx + t w - t 2 < t E J F (AT, t u t 2 , t w ) < t w (2.1) 

that characterizes the evolution of the whole system, with the time spent in the 
lower temperature partially contributing to t FFF (AT, t\, t 2 , t w ). 

Since the effects of negative cooling cycles are small, the fraction of time spent at 
T~ needs to be rather large to observe them. Fig. [I] displays a set of curves measured 
after performing negative temperature cyclings of different magnitude (thin lines) 
and compares them to four isothermal curves associated to smaller waiting times. It 
is clear from the figure that the curves associated to each temperature cycling can 
be completely superposed to a reference isothermal curve (compare, for example, 
the curve for AT = T~ — T + = -0.3K, t x = 15 min, t 2 = 1015 min and t w = 1030 
min with the isothermal curve for t w = 100 min). The system kept a memory of the 
evolution during t\ when the temperature is raised back to its original value T + . 
The time spent at the low temperature has a partial effect. There is an effective 
waiting time t FFF for the full system that verifies Eq. (|2.1| ). The effective waiting 
time t FFF decreases when one cools down the system to a lower temperature until the 
effect of the time spent at the lower temperature becomes completely negligeable. 
In the experimental curves this happens for AT = —IK. In this case, the effective 
aging time is simply i FFF = ti + t w — t 2 = 30 min. 

B. Positive temperature cycling 

In this case the straight line in Fig. [5] moves first anti-clockwise and then clock- 
wise, exposing at the end of the cycle the bold segment of the curve. 

As already discussed, the slow part of the correlation is only modified by a time- 
reparametrization and, since the temperature excursion has a short duration, it is 
hardly affected at all. 

The protocol is, from the point of view of the fast model, as follows: at time 
t = it is quenched to a subcritical temperature T~ and let age until t\. At this 
time it is taken to its critical temperature T + until time t 2 , fully reinitializing its 
age. Finally, it is quenched again into its glassy phase at time t 2 and it starts aging 
again. It is then clear that the fast model has an age t FFF = t w — t 2 < t w . 

The observed TRM curves are then a superposition of a slow part, which is 
essentially unaffected by the temperature cycle, plus a fast part that is completely 
reinitialized. The relative magnitude of fast and slow parts depends upon the height 
of the positive excursion. Thus, for given T~ , the larger T + the faster the total decay 
of the TRM as shown in Fig. |2| where the TRM curves for temperature cyclings of 
the same duration but different heights are displayed. 

Having the whole decay of the TRM curves as a function of the time- difference 
t — t w one can discuss them in great detail. One easily notices from the curves in 
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Fig. |2] that one cannot associate an effective waiting time to the full TRM decay 
after the cycling, consistently with the the additive solution obtained analytically 
for the rejuvenation process. 

Within the analytic framework, the larger the temperature T + , the smaller the 
go and C qo , implying that the fast part of the correlation and response that is 
completely refreshed (and hence the total rejuvenation) is larger, as we see in Fig. |2|. 

It has been stressed! that the TRM curve after a positive temperature cycle 
at first coincides with the isothermal TRM curve of waiting time t w ~ t — t2, the 
difference only becoming manifest after rather long times. If one compares the TRM 
curves in Fig. ^ for several values of T + with the isothermal curve for t w = 30 min, 
one observes that the departure from this reference curve is achieved later for the 
larger T + s. 

This is precisely what we would expect from the many-scales picture: Two pro- 
tocols with the same T~ and slightly different high temperatures T + and T + (as- 
sociated with g and q r ) will differ in the refreshment of the scales of C and M with 
qo < C < q' , and these are slower the smaller go- Hence, the difference in TRM 
between the two protocols will only show up at long times, the higher T + . 

This effect should become more and more marked for longer experiments. 



III. FIELD CYCLINGS 

Field cyclings at constant temperature are known to yield rather large reinitial- 
izations in the out of phase susceptibility both on increasing and on decreasing the 
Although we have not done a full analysis for the field cycling experiments, 
we may hint the reason why the arguments we used to justify a large asymmetry in 
temperature jump experiments do not carry through unaltered to the field cycling 
case. 

The effect of a time-dependent field h(t) on the dynamic equations is to add a 
term 

A(t,t') = hit) f drh(r)R(t',T) (3.1) 
Jo 



to Eq. (|1.3j) , and a term A(t,t) to Eq. (O). In the constant field situation, these 



terms have the effect of eliminating the slowest scales. The ^vsC plot of Fig. |3] for 
a constant field h coincides with the one for h = but terminates at an /z-dependent 
point (qo,Xo)- The most outstanding effect of switching on (switching off) a field is 
then to erase (create) the slowest scales between (0, x°°) an d (go, Xo)- Note that this 
is already very different from the effect of temperature changes that mainly affect 
the fastest scales in the problem. 

One could then ask if a system that has been evolving at zero field and is sud- 
dendly taken to finite h does preserve the fast parts C F and R F without reinitial- 
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ization (just the mirror image of what an increase of temperature does). In fact, 
we can propose that such a solution exists, and see where the argument takes us. 
We hence assume a separation into a fast and a slow parts of the correlation and 



response just as in Eq. (|1.7|) with the value q chosen to be the smallest possible 
correlation for a field h. 

If we now assume that on turning on the field the fast evolution remains essen- 
tially unaltered, while the slow evolution in C qo remains slow and finally tends to qo 
(as in a constant field case), we can separate out the fast part. The fast equations 
so obtained are just like Eqs. ( |1.9| )-( |iTTT| ) with additional terms which (just as in 
the temperature cycling case happened for the slow parts) are constant up to, and 
long after the field jump. The difference is that now the "bump" in this otherwise 
constant factors and terms is no longer short-lasted with respect to the fast times, 
precisely because it comes from the slow equation and from the term (|3.1| ), which 
has waiting time effects itself. 

Hence, the initial assumption that a decoupled equation for the fast part could 
be found does not hold, at least in the manner it did for the slow parts in the 
temperature cycling case. The fast parts must be altered by the field cycling in an 
important way implying no big difference between switching on and off the field. 



IV. PROPOSAL FOR NEW EXPERIMENTAL MEASUREMENTS 

The obvious direct way of constructing the x vsC plot is to perform independent 
measurements of the time-correlation of the magnetization noise C(t,t w ) and the 
integrated response x(t, t w ) at a given temperature and field. One can a posteriori 
check the independence on temperature and field of these curves. 

If the assumption of independence of temperature and field of the ^vsC curve 
in its aging part holds true, this curve can be simply obtained in a manner that 
is much more simple experimentally as it does not involve noise measurements and 
only ac susceptibilities. Indeed, let us define the quantity y(T, h) as follows 

y(T,h) = lim lim X '(",t w ) = - - • (4.1) 

LU — >U t w — >0O J 

The Edwards- Anderson parameter g s ^(T) would be then independent of the field h 
below the de Almeida-Thouless line (or "pseudo de Almeida-Thouless line" if it is 
only a dynamic crossover) and 

y(T,h) = y(T). (4.2) 
The parameter corresponds to 

q d = lim q E A(T) . (4.3) 
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Equation (|4.2| ) can be easily checked experimentally. 

Another related test is to measure the dc susceptibility 

X°°(T, h) = lim lim X (t,t w ) (4.4) 

t w — >00 t — >oo 

at constant field h, presumably obtainable from the field-cooled experimental dataj^l 
and check if it is independent of temperature everywhere below the (pseudo) 
Almeida-Thouless line. There is quite a bit of evidence in this direction in the 
literature.0 

One can then extract the curved part of the ^vsC plot of Fig. ^, that we called 
f(C) from the set of equations 

C = q d -Ty(T), 
f(C) = y(T) , (4.5) 

using T as a parameter. 

Once (and if) the validity of the approximation is verified, one can use the x vs 
C plot so obtained to check the consistency of the present reasoning for the cycling 
experiments. For example, in the positive cycling experiments the difference in 
TRM between the isothermal and the cycled experiment just after the cycle should 
be equal to the difference in height in the x vs C pl°t °f the intersection of the 
straight lines corresponding to T~ and T + (this is best seen by looking at Fig. ^). 
This non-trivial relation can be checked experimentally. 



V. CONCLUSIONS 

While it is gratifying to be able to reproduce the qualitative experimental features 
analytically, we believe one should not hasten to declare that this is evidence for 
the mean-field scenario holding strictly. Thus, infinitely separated time-scales could 
in real life be just very separated scales, the Almeida-Thouless line could be just a 
dynamic crossover, etc. 

The derivation in this paper is the first one done for a microscopic model that 
on one hand captures the main features of the experiments and on the other hand 
demonstrates the difference in the behaviour of glasses and spin-glasses. Although 
our calculation has relied on a particular model we believe that it will carry through 
to all classical mean-field spin-glass models. 

There is however a problem. A satisfying scaling of the (isothermal) experimental 
susceptibility and TRM data simultaneously was presented in Ref . []2D[] . Surprisingly 
enough, this scaling has only two time-regimes as opposed to the many time-scales 
that we here invoked to reproduce analytically the experimental results for tem- 
perature cyclings within mean-field. If the scaling used in Ref. [p9[ turns out to 
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hold for all times then the hierarchical scenario would not apply and some drastic 
modification of the present understanding must be envisaged. 
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